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Abstract 

Recently, we demonstrated rigorously the stability of charge-stripe phases in quantum- 
particle systems that are described by extended Falicov-Kimball Hamiltonians, with 
the quantum hopping particles being either spinless fermions or hardcore bosons. In 
this paper, by means of the same methods, we show that any anisotropy of nearest- 
neighbor hopping eliminates the 7r/2-rotation degeneracy of the so called dimeric and 
axial-stripe phases and orients them in the direction of a weaker hopping. Moreover, 
due to the same anisotropy the obtained phase diagrams of fermions show a tendency 
to become similar to those of hardcore bosons. 

1 Introduction 

In our recent paper [I] we addressed the problem of formation of static charge-stripe phases in 
systems of correlated quantum particles. Since an experimental evidence for the existence of 
striped phases in materials exhibiting high-temperature superconductivity has been gathered 
EI, one observes vigorous discussions of the problem of the stability of such phases in 
simple models of correlated electrons. Mainly, some Hubbard-like or t — J-like models are 
considered in this context (see a review paper by Oles |4J. Unfortunately, none of the 
methods applied for investigating these models enables one to control tiny energy differences 
between competing phases: the phase energies are calculated by means of uncontrolled 
approximations, moreover the results are biased by finite-size and boundary effects. In this 
situation, to get a deeper insight into the stability problem of striped phases, some simpler 
analogue problems |HJ El HI IHj have been studied. 

Adopting such an approach, we have proposed in Ref. pQ a kind of an extended Falicov- 
Kimball model, whose ground-state phase diagram at the half-filling and at a sufficiently 
strong coupling is amenable to rigorous investigations - a rather exceptional situation in the 
field of correlated quantum-particle systems. From a technical point of view, the standard 
spinless Falicov-Kimball model (such as that studied in [3 El) has been modified in order 
to allow for segregated phases, in the regime specified above. That is, the Hamiltonian 
has been augmented by a direct, Ising-like, nearest-neighbor (n.n.) attractive interaction 
between the immobile particles. However from the point of view of crystallization problem, 
where the immobile particles are interpreted as ions JU] > the added attractive interaction is 
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a physically motivated supplement of the hardcore repulsion, which is already incorporated 
in the spinless Falicov-Kimball model. 

We have succeeded in proving that charge-stripe phases are stable in some domains of 
phase diagrams. By varying slightly the interaction parameter of the Ising-like n.n. interac- 
tion, the system is driven from a crystalline (chessboard) phase to a segregated phase, via 
quasi-one- dimensional striped phases. 

In our studies in Ref. [T| and in this paper, we have included also a subsidiary direct 
interaction between the immobile particles, an Ising-like next nearest-neighbor (n.n.n.) in- 
teraction, much weaker than the n.n. interaction. This interaction can reinforce or frustrate 
the n.n. interaction, depending on the sign of its interaction constant. We find it useful in 
discussing similarities and differences between the cases of spinless fermions and hardcore 
bosons. 

The obtained results on the stability of striped phases open new possibilities of investigat- 
ing various characteristics of these phases. Since, it is the role of a hopping anisotropy that 
has received a considerable attention in the recent literature of the subject jU HH [T2J 1TB"] . 
we have turned to this problem, and in this paper we discuss the role of a weak anisotropy 
of n.n. hopping. 

The paper is organized as follows. In next section we define the model studied in this 
paper, describe symmetries of its grand-canonical phase diagrams, and introduce strong- 
coupling expansions of the ground-state energy (effective Hamiltonians) and the correspond- 
ing phase diagrams. Then, in Section 3, we carry out a detailed analysis of phase diagrams 
due to truncated effective Hamiltonians of orders not exceeding four, in presence of a weak 
anisotropy of n.n. hopping. In Section 4, we draw conclusions and provide a summary. 
Various technical details a placed in Appendices A, B, and C. 

2 The model and its ground-state phase diagram — 
basic properties 

2.1 The model 

The model to be studied is a simplified version of the one band, spin 1/2 Hubbard model, 
known as the static approximation (one sort of electrons hops while the other sort is immo- 
bile), augmented by a direct Ising-like interaction V between the immobile particles. Only 
a n.n. hopping is taken into account and we allow for its anisotropy. The total Hamiltonian 
of the system reads: 



H = H FK + V, (1) 

h fk = -t h ( c * c y + c v Cx ) ~ tv ( c * c y + c v Cx ) +U Y1 \ c * Cx ~ \) Sx > ( 2 ) 

{x,y)i,h (&>!/) i,« x 

In the above formulae, the underlying lattice is a square lattice, denoted A, consisting of 
sites x, y, . . ., whose number is |A|, having the shape of a a/|A[ X ypM torus. In (|2I3|) and 
below, the sums Y^(xy) ih ' or ^(a; !/)»«' * = 1>2>3' stand for the summation over all the i-th 
order n.n. pairs of lattice sites in A, directed horizontally (h) or vertically (v), with each 
pair counted once. 
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The subsystem of mobile spinless particles (here-after called the electrons) is described 
in terms of creation and annihilation operators of an electron at site x, c+, c x , respectively, 
satisfying the canonical anticommutation relations (spinless fermions) or the commutation 
relations of spin 1/2 operators S£,S~,S* (hardcore bosons). The total electron-number 
operator is N e = J2 x c t c x, and (with a little abuse of notation) the corresponding electron 
density is p e = iV e /|A|. There is no direct interaction between the mobile particles. Their 
energy is due to the anisotropic hopping, with t^ (t v ) being the n.n. horizontal (vertical) 
hopping intensity, and due to the interaction with the localized particles, whose strength is 
controlled by the coupling constant U . 

The Hamiltonian Hfk is well-known as the Hamiltonian of the spinless Falicov-Kimball 
model, a simplified version of the Hamiltonian put forward in [Hj. A review of rigorous 
results concerning this model and an extensive list of relevant references can be found in 

psjQuini). 

The subsystem of localized particles (here-after called the ions), is described by a col- 
lection of pseudo-spins {s x } x€A , with s x = 1, —1 (s x = 1 if the site x is occupied by an 
ion and s x = —1 if it is empty), called the ion configurations. The total number of ions 
is Ni = ^2 x (s x + l)/2 and the ion density is pi = A^/|A|. In our model the ions interact 
directly by means of the isotropic Ising-like interaction V. 

Clearly, in the composite system, whose Hamiltonian is given by with an arbitrary 
electron-ion coupling U, the particle-number operators N e , N iy and pseudo-spins s x , are 
conserved. Therefore, the description of the classical subsystem in terms of the ion config- 
urations S = {s x } xeA remains valid. Whenever periodic configurations of pseudo-spins are 
considered, it is assumed that A is sufficiently large, so that it accommodates an integer 
number of elementary cells. 

2.2 The ground-state phase diagram in the grand-canonical en- 
semble 

In what follows, we shall study the ground-state phase diagram of the system defined by (JTJ) 
in the grand-canonical ensemble. That is, let 

H (p e , p^ = H - p e N e - p,iNi, (4) 

where p e , pi are the chemical potentials of the electrons and ions, respectively, and let 
Es{pe,Pi) be the ground-state energy of H(p e ,pi), for a given configuration S of the 
ions. Then, the ground-state energy of H(p e ,pi), Eg (p e , Pi), is defined as Ea(p e ,Pi) = 
mm {Es (p e , Pi) '■ S}. The minimum is attained at the set G of the ground-state configu- 
rations of ions. We shall determine the subsets of the space of the Hamiltonian's energy 
parameters, where G consists of periodic configurations of ions, uniformly in the size of the 
underlying square lattice. 

In studies of grand-canonical phase diagrams an important role is played by unitary 
transformations (hole-particle transformations) that exchange particles and holes, c x c x — > 
1 — c x c x , s x — > —s x , and for some (p° e1 p®) leave the Hamiltonian H (p e , p^ invariant. For 
the mobile particles such a role is played by the transformations: c+ — > e x c x , where e x = 1 
for bosons while for fermions e x = 1 at the even sublattice of A and e x = — 1 at the odd 
one. Clearly, since Hq is invariant under the joint hole-particle transformation of mobile 
and localized particles, H (p e , p^ is hole-particle invariant at the point (0, 0). At the hole- 
particle symmetry point, the system under consideration has very special properties, which 
simplify studies of its phase diagram Moreover, by means of the defined above hole- 
particle transformations one can determine a number of symmetries of the grand-canonical 
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phase diagram [18]. The peculiarity of the model is that the case of attraction (U < 0) 
and the case of repulsion (U > 0) are related by a unitary transformation (the hole-particle 
transformation for ions): if S = {s x } xeA is a ground-state configuration at (fa, fa) for U > 0, 
then — S = {— s x } xeA is the ground-state configuration at (/i e , —fa) for U < 0. Consequently, 
without any loss of generality one can fix the sign of the coupling constant U. Moreover 
(with the sign of U fixed), there is an inversion symmetry of the grand-canonical phase 
diagram, that is, if S is a ground-state configuration at (/i e , fa), then — S is the ground-state 
configuration at (— fa,— fa). Therefore, it is enough to determine the phase diagram in a 
half-plane specified by fixing the sign of one of the chemical potentials. 

Our aim in this paper is to investigate the ground-state phase diagrams of our systems, 
for general values of the energy parameters that appear in H (/i e , fa). According to the state 
of art, this is feasible only in the strong- coupling regime, i.e. when \t/U\ is sufficiently small. 
Therefore, from now on we shall consider exclusively the case of a large positive coupling U, 
and we express all the parameters of if (/z e ,/Zj) in the units of U, preserving the previous 
notation. 

2.3 The ground-state energy and phase diagram in the strong- 
coupling regime 

In the strong-coupling regime, the ground-state energy E$ (/i e , fa) can be expanded in a 
power series in t. One of the ways to achieve this, for fermions and for hardcore bosons, is 
a method of unitarily equivalent interactions [TP] . The result, with the expansion terms up 
to order four shown explicitly, reads: 
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in the case of hopping fermions, and 

E b s (ve,Vi) = {E b s f ] (a*., A*) + {R b sf\ 
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in the case of hopping hardcore bosons. The both of the above expressions are given up to 
a term independent of the ion configurations and the chemical potentials. In (JHJ and (jHJ), 
we have set th = t, and t v = y/jt, with < 7 < 1, P denotes a (2 x 2)-sites plaquette of 
the square lattice A, sp stands for the product of pseudo-spins assigned to the corners of 

P. The remainders, ( R s J and [Rs) , which are independent of the chemical potentials 

and the parameters W and e, but dependent on 7, collect those terms of the expansions 
that are proportional to t 2m , with m = 3, 4, . . .. It can be proved that the above expansions 

< 1 - 16t 



are absolutely convergent, uniformly in A, provided that t < 1/16 and \ji e 
Moreover, under these conditions the ground-state densities of particles satisfy the half-filling 
relation: p e + pi — 1. 

We note that on taking into account the inversion symmetry of the phase diagram in the 
(fi e , fii) -plane and the fact that the ground-state energies depend only on the difference of 
the chemical potentials, in order to determine the phase diagram in the stripe \/j, e \ < 1 — 16t 



it is enough to consider the phase diagram at the half-line \x e 
this half-line we set //j = \i. 



0, fa < (or fii > 0). At 



3 Phase diagrams according to truncated expansions 

Due to the convergence of the expansions © and (jOJ), it is possible to establish rigorously 
a part of the phase diagram, that is the ground-state configurations of ions are determined 
in the space of the Hamiltonian's energy parameters everywhere, except some small regions. 
This is accomplished by determining the ground-state phase diagram of the expansion trun- 
cated at the order k {the k-th order phase diagram), that is according to the fc-th order 
effective Hamiltonians (Eg)^ (0, p) and {E h s )^ (0, p). The remainder of the expansion, 
that consists of a term of the next order, k' > k, and other higher-order terms, cannot 
modify the obtained phase diagram, except some narrow regions along the phase boundaries 
whose width is of the order k! . The procedure of constructing the phase diagram is recursive: 
the phase diagram of the effective interaction truncated at the order k' is constructed on the 
basis of the phase diagram obtained at the preceding order k: the conditions imposed on the 
ground-state configurations by the k-th. order effective Hamiltonian have to be obeyed by 
the ground-state configurations of the k'-th order effective Hamiltonian. In other words, the 
k'-th order terms of the expansion cannot change the hierarchy of configuration's energies 
established by the k-th order effective Hamiltonian; they can only split the energies of con- 
figurations in cases of degeneracy. At each step of the construction, we use the m-potential 
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method introduced in [2JJ, with technical developments described in [TBI 1221 120] • For details 
specific to this paper see also Ref. [JJ and Appendix A. 

In Ref. [JJ we have obtained the ground-state phase diagrams, according to the fourth- 
order isotropic effective Hamiltonians (Eg)^ (0, fi) and (E s )^ (0, fi) | , for hole-particle 

7=1 ^~ 

symmetric systems (/i = 0) with a weak (of fourth order) n.n.n. subsidiary interaction, and 
for unsymmetrical systems without the subsidiary interaction (e = 0), see the top phase 
diagrams in Figs. I1J2J31 Second-order phase diagrams of Ref. [JJ consist exclusively of phases 
whose configurations are invariant with respect to 7r/2-rotations, with macroscopic degen- 
eracies (i.e. the number of configurations grows exponentially with the number of sites) 
at the boundaries of phase domains. Stripe phases, whose configurations are not invariant 
with respect to 7r/2- rotations, appear on perturbing the second-order phase diagrams by the 
fourth-order isotropic interactions. Here, we would like to observe the influence of a weak 
anisotropy of n.n. hopping on these stripe phases. Since we are working with truncated 
effective Hamiltonians, we have to assign an order to the deviation of the anisotropy param- 
eter 7 from the value 1 (corresponding to the isotropic case). Therefore, we introduce an 
anisotropy order, a, and a new anisotropy parameter, /3 a : 

7 = 1- dat a . (7) 

The orders of the deviations can be neither too small, not to modify the second-order effective 
Hamiltonians, nor too high, to effect the considered effective Hamiltonian of the highest 
order. Since here, the highest order of the effective Hamiltonians is k = 4, the weakest 
admissible deviation from the isotropic case corresponds to the highest anisotropy order a = 
2. Then, we can consider an intermediate deviation, i.e. < a < 2. The strongest deviation, 
i.e. a = 0, is not admissible, since it modifies the second-order effective Hamiltonians. 

3.1 The smallest deviation from the isotropic case 

In the sequel, we drop the arguments of ground-state energies, that is we set (j^gj (0, /i) = 

/ ,\(4) 

( E J S ) , etc. In the case of the smallest deviation from the isotropic case the fourth-order 
effective Hamiltonian for fermions reads: 



A (4) = f R /V 4) 



H = [E s 



) 



t 4 

#2-7 22 SxS ^ ( 8 ) 



4 

7=1 (x,vh,v 



while for hardcore bosons only the first term, representing the isotropic fourth-order effective 
Hamiltonian, has to be changed properly. 

Apparently, the effective Hamiltonians of order zero and two are isotropic, and there is 
no difference between the cases of hopping fermions and hopping bosons. The corresponding 
phase diagrams, determined in Ref.[JJ, are built of three phases, i.e. sets of configurations 
with the same energy in some open regions of phase diagrams, called phase domains, S + , 
S-, and S c b- The first two phases consist of single translation invariant configurations, 
one completely filled with ions and one completely empty, respectively; a representative 
configuration of the chessboard phase, S c b, is shown in Fig. The three named phases 
coexist at the open half-line W = —2t 2 , fi — 0, e > 0. However, if e attains zero, the finite 
(independent of the size of A) degeneracy of the open half-line changes into a macroscopic one. 
It is the point W = —2t 2 , fi = 0, e = around which the fourth-order diagrams, containing 
various stripe phases, have been constructed in Ref. 0. Now, to observe the effect of the 
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fourth-order anisotropy term of (JHJ) the construction of the fourth-order diagrams has to be 
carried out again. As in the isotropic case, this is facilitated by introducing new variables, 
uj, 6, and e: 

W = -2t 2 + t 4 cu, p = t 4 5, e = t 4 e, (9) 
and rewriting the fourth order effective Hamiltonian in the form, 
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with analogous expressions in the bosonic case, and with the isotropic potentials ( W T \ 



7=1 



given by 



(4) 



7=1 



= -5 (s 5 + 1) + — I w - - J + — (3 - e) + 



7=1 



(12) 



b\(4) 



7=1 



= -5 (s 5 + 1) + — (u - - J s * s y + 32 ( 5 ~ £ ) S s ^ + 

1 vW' 1 vW' . . 



<^,2/>3 



(13) 



In the above formulae T stands for a (3 x 3)-sites plaquette of a square lattice, later on 
called the T '-plaquette, and s$ stands for the central site of a T-plaquette. The summation 
in double-primed sums runs over a T-plaquette. 

We would like to get an idea of the phase diagram in the space of the four energy 
parameters (u, e, 02, that appear in the Hamiltonian. Due to the fact that the domains 
occupied by the phases are polyhedral sets (see Appendix A), this goal can be achieved by 
studying phase diagrams in two-dimensional hyperplanes. Of particular interest are those 
hyperplanes that result from intersecting the four- dimensional space by hyperplanes 5 = 
and p 2 — const, and by hyperplanes e = and /3 2 = const. A collection of such sections for 
suitable values of /?2, forming a finite increasing from zero sequence, enable us to observe how 
the anisotropy effects our system if it is hole-particle symmetric and if it is not, respectively. 

Specifically, in Fig. ^ we show phase diagrams in the hole-particle-symmetric case, while 
in the absence of the hole-particle symmetry, the phase diagrams are shown in Fig. El - for 
fermions, and in Fig. |3] for bosons. 

The phases that appear in these phase diagrams can conveniently be described in terms 
of the phases found in ground-state phase diagrams of the isotropic systems pQ, with Hamil- 
tonian H \ v and presented in Fig. H21 The phases of anisotropic systems, considered here, 
either coincide with or are simple modifications of the isotropic phases. In fact, in Fig.lT2lwe 
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display only representative configurations of the phases of isotropic systems. The remaining 
configurations can be obtained by means of the spatial symmetry operations of Hq\ x , like 
translations and rotations by tt/2. The numbers in curly brackets, placed by the symbols 
of phases, stand for the numbers of the T-plaquette configurations (according to Fig. ITlj) 
that are obtained by restricting the configurations of a phase to a T-plaquette. By the same 
symbols as the phases we denote also their domains. 

Among the phases of isotropic systems, we can distinguish a class of dimeric phases, Sdi, 
. . . , <Sd4, and a class of axial-stripe phases, S v \, S v2 , and S v3 . The configurations of dimeric 
phases consist of isolated pairs of n.n. occupied sites (in the sequel called the dimers). In 
the configurations of axial-stripe phases, the ions fill completely some, parallel to one of 
the axes, lattice lines, so that a periodic pattern of stripes is formed. Out of the dimeric or 
axial-stripe configurations of Fig. ^1 only those with dimers or stripes, respectively, oriented 
vertically appear in the phase diagrams of anisotropic systems. Such a restricted phases are 
marked in anisotropic phase diagrams by the additional superscript, v. 

All the phases presented in Fig. H21 except Sd2 and SdA, contain exclusively periodic con- 
figurations, related by the spatial symmetries of H \^ =1 , hence of the degeneracy independent 
of the size of the underlying lattice. The set S pc b (of plaquette-chessboard configurations) 
contains configurations built of elementary plaquettes with occupied sites, forming a square 
lattice according to the same rules as filled sites form a square lattice in the chessboard 
configurations of S c b- The remaining phases consist of configurations that have a quasi-one- 
dimensional structure. That is, they are built of sequences of completely filled with ions 
lattice lines of given slope. Such a configuration can be specified by giving the slope of the 
filled lattice lines of a representative configuration and than the succession of filled (f) and 
empty (e) consecutive lattice lines in a period. For instance, the representative configuration 
of Sdi (see Fig. H2J) is built of filled lines with the slope 2 and, in the period, two consecu- 
tive filled lines are followed by two empty lines, which is denoted (2;2/, 2e). This kind of 
description of the remaining quasi-one-dimensional configurations is given in Fig. The 
phase Sdd is an example of a diagonal- stripe phase. 

Only in the domains Sd2 and SdA, which appear in the phase diagrams shown in Figs. El 
that is off the hole-particle symmetry case, the situation is more complex, their degeneracy 
grows indefinitely with the size of the lattice. In iS rf2 one can distinguish two classes, Sd2a 
and Sd2b-, of periodic configurations with parallelogram elementary cells. A configuration 
in Sd2a consists of vertical (horizontal) dimers of occupied sites that form a square lattice, 
where the sides of the elementary squares have the length 2\/2 and the slope ±1. In a 
configuration of Sd2bi the elementary parallelograms formed by dimers have the sides of the 
length 2\/2 and the slope ±1, and the sides of the length vTO and the slope ±1/3. Two 
configurations, one from Sd2a and one from Sd2b, having the same kind of dimers (vertical 
or horizontal), can be merged together along a "defect line" of the slope ±1 (the dashed 
line) , as shown in Fig. without increasing the energy. By introducing more defect lines 
one can construct many ground-state configurations whose number scales with the size of 
the lattice as exp {constVA). The same kind of degeneracy is in SdA- This phase consists 
of three classes of periodic configurations of dimers. In the class SdAa, the elementary cell 
can be chosen as a parallelogram whose two sides of length 3 are parallel to dimers (which 
are vertical or horizontal). If the dimers are oriented vertically, then the other two sides 
have the slope 1/2 and the length a/5. By reflecting an elementary cell of SdAa in a lattice 
line passing through its side that is parallel to dimers, we obtain an elementary cell of the 
class SdAb- In the third class, SdAc, an elementary cell can be chosen as a rhomb formed 
by the centers of dimers, with the sides of length vTO. Two configurations, one from Sd& a 
and one from SdAb, having the same kind of dimers (vertical or horizontal), can be merged 
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together along a " defect line" parallel to dimers (dashed line in Fig. IT3j) without increasing 
the energy. In this way, a numerous family of configurations can be constructed, with the 
number of configurations growing like exp (consty/X). 

It follows from the polyhedral shape of the phase domains, that the set of values of fa 
is partitioned into open intervals, where the boundary lines of phase domains do not change 
their direction, only their distance to the origin varies in an affine way. The boundary points 
of these open intervals are special values of the anisotropy for the phase diagrams, and in the 
sequel we call them the critical values. As a critical value of anisotropy is approached, some 
boundary lines merge into a line or a point, which results in disappearance of some phase 
domains. And vice versa, some points and lines break off, creating new phase domains. 

In particular, in the hole-particle symmetric case it can be inferred from (Fig. that at 
least up to fa = 10 - for hopping fermions, and at least up to fa = 5 - for hardcore bosons, 
there is only one critical value of fa. For fermions, it amounts to fa = 7, where the phase 
Sdd disappears, while for hardcore bosons it is fa = 2, where the phase S pc b disappears. 

Off the hole-particle symmetry, for bosons, there are no critical values of fa, at least up 
to fa — 3. On the other hand, for fermions and for fa < 6 there are as many as three critical 
values of fa. In increasing order, the first is fa — 1, where the phases S% x , S^, and S% 2 
appear. The second is fa = 2, where the phases S^ 3 , S^, and S^ 3 appear. And the last one 
is fa = 3, where Sdd is replaced by S% 2 - 

If there are no more critical values of the anisotropy parameter fa, the shape of the fourth- 
order phase diagrams for any value fa larger than the greatest considered in Figs, d EH 
remains the same. On increasing the anisotropy parameter, the diagrams undergo only 
some translations. To verify whether for the values of fa larger than those considered in this 
section new critical values do appear, we proceed to investigating stronger, i.e. intermediate, 
deviations from the isotropic case. 

3.2 The intermediate deviation from the isotropic case 

Now, the fourth-order effective Hamiltonian for fermions assumes the form: 



t 2+a 



4 

7=1 (*,»>!,« 



and for hardcore bosons an analogous formula holds true. Obviously, the phase diagrams 
in the zeroth and second orders remain the same as in the isotropic case, described above. 
Therefore, we proceed to constructing the phase diagram in next order, which is (2 + a)-order 
with < a < 2, and the corresponding effective Hamiltonian reads: 

Ef^=Ef-p^ £ SxSy , (15) 

(2) 

where Eg stands for the, common for fermions and bosons, second-order effective Hamil- 
tonian. For the reasons given in the previous subsection, we consider a neighborhood of the 
point W = —2t 2 , /i = 0, £ = 0, where the energies of all the configurations are equal Ref.pQ. 
In this neighborhood it is convenient to introduce new variables, 5' , e', and a/, 



t 2+a 5', e = t 2+a e', W = -2t 2 + t 2+a u', (16) 
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Figure 1: The case of the weakest anisotropy and of hole-particle symmetric systems (/x = 0). Phase diagrams of (h£\ 

(given by 1111 ') — left column, and (Hip) — right column, for an increasing sequence of values of 0%. The representative ion 
configurations of the displayed phases are shown in Fig, 1121 (for more comments see Section 3). For fermions, the critical value 
is P2 = 7, while for bosons it is P2 = 2. The equations defining the boundary lines of the phase domains are given in Tab. HI 
and Tab. l2*lof Appendix B, while the corresponding zero-potential coefficients {«;} — in Tab. 181- Tab. U5l of Appendix C. In 
the bottom diagrams, the regions surrounded by dotted circles, are reconsidered in the case of an intermediate anisotropy. 
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Figure 2: The case of the weakest anisotropy, off the hole-particle symmetry, with e = 0. The phase diagram of (b£\ 

(given by llll l for an increasing sequence of values of 02- The representative ion configurations of the displayed phases 
are shown in Fig. 1121 (for more comments see Section 3). The critical values are: 02 = 1, 02 = 2, and 02 = 3. In the 
blank region of the phase diagram for the critical 02 = 1, the following T-plaquette configurations have the minimal energy: 
26,52,53,68,81,88,101,116,117,143 (see Fig. 1141 . By means of these T-plaquette configurations one can construct SJJ 2 i ^ddi an d 
many other configurations. The equations defining the boundary lines of the phase domains are given in Tab. l3l of Appendix 
B, while the corresponding zero-potential coefficients {cti} — in Tab. 1161 — Tab. !22^ of Appendix C. In the bottom diagram, the 
regions surrounded by dotted ellipses are reconsidered in the case of an intermediate anisotropy. 
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Figure 3: The case of the weakest anisotropy, off the hole-particle symmetry, with e = 0. The phase diagram of (H^,) 
(given by llll ); the isotropic diagram (ySg = 0) and the anisotropic diagram (@2 = 3). The representative ion configurations 
of the displayed phases are shown in Fig. 1121 (for more comments see Section 3). Here, no critical values have been detected. 
The equations defining the boundary lines of the phase domains are given in Tab. |^ of Appendix B, while the corresponding 
zero-potential coefficients {oj;} — in Tab. I2ijl and Tab. 1241 of Appendix C. In the right diagram, the regions surrounded by 
dotted ellipses are reconsidered in the case of an intermediate anisotropy. 

and rewrite the expansion up to the order 2 + a: 



E( s +a) = l — \ ~ 5 ' ( s * + x ) + T SxS y SxS v ~ y 2 




S X Sy 



{x,y)2 (x,y) 



1,1) 



(17) 



where 



Hp ^ — ^ (S X + 1) + ^ S X Sy ^ S X Sy ^ ^ S X Sy. 

(as > y)i (x,v)2 (x,y)i,v 



The summations in the primed sums run over a plaquette P. The plaquette potentials Hp*' 
have to be minimized over all the plaquette configurations. As in the previous case, we are 
interested in phase diagrams of hole-particle symmetric systems (5' = 0) or unsymmetrical 
systems with e' — 0. It turns out that for such energy parameters and plaquette configu- 
rations the potentials Hp +a ^ are m-potentials. The resulting phase diagrams are shown in 
Fig. Eland Fig. We note that each of the points a: u' = f3 a , e' = f3 a , b: J = 0, 5' = 0, 
c~: to' — 2j3 a , 5' = —ft a , and c + : to' — 2(3 a , 5' = j3 a , is the coexistence point of three periodic 
phases. That is, the only plaquette configurations minimizing the potential Hp +a ^ at such a 
point are those obtained by restricting the configurations of coexisting phases to a plaquette 
(these plaquette configurations are shown in Fig. 0| 

Now, following our recursive procedure of constructing phase diagrams to some order, we 
are ready to investigate the effect of fourth-order interactions. As in earlier steps, it is enough 
to consider neighborhoods of the coexistence points a, b, and c~ of Figs. 14151 (the diagram 
in a neighborhood of c + can be obtained from that around c~ by a symmetry operation). 

A. A neighborhood of a 

Here 5' = 0, hence the systems are hole-particle invariant. A convenient change of variables 
is: 

J = (3 a + t 2 - a to, e' = p a + t 2 ~ a e. (19) 
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Figure 4: The case of an intermediate anisotropy (0 < 
a < 2) and of hole-particle symmetric systems (/i = 0). The 
phase diagram (common for hopping fermion and hardcore 
boson systems) of ffp + °'. The coordinates of point a are 
lu' = f3 a , e' = Pa- The boundary lines of iS^ 2 > from left 
to right, are: e' = u>' and e' = — u/ + 2j3 a - The boundary 
between iS se9 and S co is u>' = j3 a - 



Figure 5: The case of an intermediate anisotropy (0 < 
a < 2), off the hole-particle symmetry, with e' = 0. The 
phase diagram (common for hopping fermion and hardcore 
boson systems) of Hp +a \ In the (id', <5')-plane, b = (0,0), 
c+ = (2/3 a ,/3 a ), c~ = (2j8 ,-/3 ). The boundary lines of 
<S- , from left to right, are: 5' = 0, 8' = — ^- and <5' = —a/ + 
/3 a . The boundary lines of 5+ are obtained by changing 
5' -» -5'. 



In these variables, the fourth-order effective Hamiltonian reads: 



(i) 



t 2+a 
IT 



(2+o) 

p 



5'=0 

£'=<d'=/3 a 



with a similar formula for hard-core bosons, where 



4) 



(i) 



7=1 



(20) 



H 



(2+o) 



A* 



<5'=0 
e'=o/=/3 a 



S X Sy 



(21) 



In the previous order, it has been established that the only plaquette configurations min- 

l P 



■ ■ ■ Tj(2+a) 

lmizmg H 



y_ are ^ e ones obtained by restricting to a plaquette P the periodic 

e'=J=Pa 

configurations S + , iS>_, S c b, and S% 2 . Let us denote this set of plaquette configurations by 



Sp. Consequently, the minimization of the fourth-order potentials ( H^, 



(i) 



> ( H t) 



7=1 



7=1 



should be carried out only over the set S T of those T-plaquette configurations whose restric- 
tion to a plaquette P belongs to Sp. The configurations of the set S T are displayed in Fig. Hoi 



It appears that on the set <Sf>, the potentials (h^ 



(4) 



7=1 



7=1 



are m-potentials. 



The obtained phase diagrams, independent of the anisotropy parameter (3 a , are shown in 
Fig. El and Fig. □ 



B. A neighborhood of b 

Here the system is not hole-particle symmetric and e' = 0. A convenient change of variables 
is: 



J = t 2 ~ a u, 5' = t 2 ~ a 5. 



(22) 
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Figure 6: The case of an intermediate anisotropy (0 < 
a < 2) and of hole-particle symmetric systems (fi = 0). 
The phase diagram of I E J S I in a neighborhood of point 
a (u)' = f3 a , e' = f3a) (see formula I2UI V The equations 
defining the boundary lines of the phase domains are given 
in Tab, 151 of Appendix B. 



s 



10 -i 




Figure 7: The case of an intermediate anisotropy (0 < 
a < 2) and of hole-particle symmetric systems (fi = 0). 
The phase diagram of (-E5) in a neighborhood of point 
a (li/ = f3 a , e' = f3 a ) (see formula 1201 '). The equations 
defining the boundary lines of the phase domains are given 
in Tab. 151 of Appendix B. 



Then, the fourth-order effective Hamiltonian takes the form, 



(4) 



(2+a) 
P 



<5'=0 
s'=w'=0 



~2 



J 



7=1 



where 



H 



(2+a) 



<5'=0 
e'=u>'=0 



4 ^ 

(x,y)i,v 



S X Sy 



with a similar formula for hard-core bosons. The minimum of H 



(2+a) 



8'=0 
e'=uj'=0 



(23) 



(24) 



is attained at 



the configurations belonging to Sp, i.e. the plaquette configurations obtained by restricting 
the periodic configurations S + , <S_, and S^ 2 to a plaquette P. Let be the corresponding 
set of T-plaquette configurations (there are no vertical pairs of n.n. sites occupied by one 

(4) 

are not the m-potentials. 



ion). This set is shown in Fig.^J Here the potentials \ H, 
The obtained phase diagrams are shown in Figs. |Hl El 



7=1 



C. A neighborhood of c 

Here the system is not hole-particle symmetric and e' = 0. A convenient change of variables 
is: 



u' = 2(3 a + t 2 ~ V 5' = -p a + t 2 - a 5. 
The fourth-order effective Hamiltonian reads: 

(4) t 2+a 



E f 



E« 



(2+a) 
P 



4eK 

S' = -l3 a ,L0' = 2l3 a * T 



(4) 



(25) 



(26) 



7=1 
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with a similar formula for hard-core bosons, where 



(2+0) 



s'=0 = X + % + + !) ( 27 ) 

6'=-p a ,u'=2p a (x,y)i, h 



The potentials if 



(2+a) 



£ ,_ are minimized by restrictions to a plaquette P of periodic 

6'=-l3 a ~u'=2l3 a 



configurations 6L, <S C &, and S% 2 , that constitute the set Sp . The corresponding set of 
T-plaquette configurations consists of configurations where no 

/ A (4) 

occupied by two ions (see Fig.lTTj). Here the potentials I Hj,\ 

7=1 

The corresponding phase diagrams are shown in Fig. and Fig. ^2 



rorizontal pair of n.n. sites is 
are not the m-potentials. 



4 Discussion of the phase diagrams and conclusions 

In subsection 3.2, we have obtained phase diagrams according to the truncated effective 
Hamiltonians of fourth order, with an intermediate anisotropy of hopping. These diagrams 
have been constructed to see what happens to the fourth-order phase diagrams with the 
weakest anisotropy if the anisotropy parameter j3 2 grows beyond the values taken into ac- 
count in subsection 3.1. Let us consider firstly the diagrams in a neighborhood of point a. 
Apparently, up to a rescaling the phase diagram in Fig. |U] looks the same as the part in the 
dotted circle of the fermionic phase diagram in Fig. ^ Similarly, up to a rescaling the phase 
diagram in Fig. [7| looks the same as the the part in the dotted circle of the bosonic phase 
diagram in Fig. ^ The same remarks apply to the phase diagrams in neighborhoods of points 
b and c~. The phase diagram in Fig. |S] reproduces the part in the upper dotted ellipse in 
Fig. |21 while the diagram in Fig. El- the part in the upper dotted ellipse in Fig. El Then, the 
phase diagram in Fig. El reproduces the part in the lower dotted ellipse in Fig. |21 while the 
diagram in Fig. El the part in the lower dotted ellipse in Fig. El Therefore, we conclude 
that in the case of the weakest anisotropy in the fourth-order effective Hamiltonians, we 
have determined all the critical values of the anisotropy parameter (3 2 . For the both kinds 
of hopping particles, if (3 2 exceeds the greatest critical value, the obtained phase diagram 
undergoes only a translation (varying with f3 2 ). 

Now, the basic question to be answered is concerned with the relation between these 
fourth-order phase diagrams and the phase diagrams of quantum systems described by 
Hamiltonians H Q . By adapting the arguments presented in Refs. [221120]; we can demon- 

rrm ( f h \ (4) 

strate, see for instance Ref. j23|, that if the remainders, ( R s J and \Rg) , are taken into 
account, then there is a sufficiently small to such that for t < t the phase diagrams of quan- 
tum systems look the same as the phase diagrams according to the effective Hamiltonians 
truncated at the fourth order, except some narrow regions, of width 0(t 2 ) (at the diagrams 
displayed above), located along the phase-domains boundaries, and except the domains S^ 2 
and S% 4 . For t < to and for each domain Sz>, which is different from S^ 2 and <S^ 4 , there 
is a nonempty two-dimensional open domain that is contained in the domain So and 
such that in Sp the set of ground-state configurations coincides with Sp. Moreover, in com- 
parison with the critical values of the anisotropy parameter j3 2 , determined according to the 
fourth-order effective Hamiltonians, the corresponding critical values of the quantum systems 
described by Hamiltonians H Q differ by 0(t 2 ), i.e. for a quantum system 7 = 1 — (3 2 t 2 + 0{t A ). 

Remarkably, in the fourth-order the hole-particle symmetric phase diagrams of fermions 
and of hardcore bosons are geometrically similar. That is, a phase diagram of hard core 
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Fi gUre 8; The case of an intermediate anisotropy 

(0 < a < 2), off the hole- particle symmetry, with s' = 0. 
/ f \(4) 

The phase diagram of I Ei, I in a neighborhood of point 
b (lu' = 0, S' = 0), see formula 1231 , The equations defining 
the boundary lines of the phase domains are given in Tab.lol 
of Appendix B, while the corresponding zero-potential co- 
efficients {oti} in Tab, !25l of Appendix C. 
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Figure 10: The case of an intermediate anisotropy 

(0 < a < 2), off the hole- particle symmetry, with e' = 0. 

„ / , \ (4) 

The phase diagram of I E s j in a neighborhood of point 
c (u)' = 2f3 a , S' = —fia), see formula 1261 , The equations 
defining the boundary lines of the phase domains are given 
in Tab. |7| of Appendix B, while the corresponding zero- 
potential coefficients {«;} in Tab. !2Yl of Appendix C. 



-6 a -4 -2 b 2 4 




8 



Figure 9: The case of an intermediate anisotropy 
(0 < a < 2), off the hole- particle symmetry, with e' = 0. 
The phase diagram of (Eg) ^ in a neighborhood of point b 
(a/ = 0, 5' = 0)), see formula 1231 , The equations defining 
the boundary lines of the phase domains are given in Tab.lol 
of Appendix B, while the corresponding zero-potential co- 
efficients {cti} in Tab. l2(jl of Appendix C. 



8 



4 - 




Figure 11: The case of an intermediate anisotropy 
(0 < a < 2), off the hole- particle symmetry, with e' = 0. 
The phase diagram of (Sg)' 4 ' in a neighborhood of point 
c~ (lu 1 = 2{3 a , 8' = — j3 a ), see formula 1261 . The equations 
defining the boundary lines of the phase domains are given 
in Tab. |7| of Appendix B, while the corresponding zero- 
potential coefficients {cti} in Tab. !28l of Appendix C. 
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Figure 12 ! The representative configurations of the phases that appear in the phase diagrams of Ho\^_ 1 . The remaining 
configurations can be obtained by applying the spatial symmetries of Hq\ As a representative configuration of the set 5^2) 

whose degeneracy grows like exp (const\/~K) , we show a configuration with one defect line (the dashed line) ; the continuous line 
is a guide for the eye. Explanations of the symbols on the top of the representative configurations and other comments are 
given in Section 3. 
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S M (prl/3) 

Figure 13: The representative configurations of the set whose degeneracy grows like exp (const\/A). The left configu- 
ration is an example of configurations with defect lines. Here, there is one defect line (the dashed line, the continuous line is a 
guide for the eye): the vertical lattice line separating a periodic configuration of vertical dimers from its vertical translate by 
one lattice constant. The right configuration is a periodic configuration of dimers. For more comments see Section 3. 



bosons, with any fa > 0, can be obtained from a phase diagram of fermions, with fa — 5, 
by the translation whose vector reads: uj = —7, e = —3, and fa = —5. The existence of this 
translation vector is related to the fact that for both kinds of systems there is one critical 
value of the anisotropy parameter. Additionally, for fermions with fa < 7, it is necessary 
to replace the phase Sdd in the central domain by S pc i,. Off the hole-particle symmetry, the 
relation between the bosonic and fermionic phase diagrams is not that close. For fermions, 
there are three critical values of fa, while for hardcore bosons there is no critical values. Thus, 
the system with hopping hardcore bosons is less sensitive to the anisotropy of hopping, than 
the system with hopping fermions. Nevertheless, if e = 0, then the phase diagrams of both 
kinds of systems are topologically similar, except that in the bosonic phase diagrams the 
phases S% 3 and S 2 are missing. However, we know from Ref. [22] that this deficiency can be 
removed by switching on the n.n.n. interactions with negative e. 

In the fourth-order effective Hamiltonians (jHJ), the weakest anisotropy of n.n. hopping 
assumes the form of a fourth-order attractive n.n. interaction in vertical direction (i.e. the 
direction of a weaker hopping). This interaction favors n.n pairs of occupied or empty 
sites that are oriented vertically. As a result, the dimeric and axial-stripe phases oriented 
vertically are stabilized for any value of the anisotropy parameter fa, while S pc b and Sdd 
are replaced by S% 2 above a critical value of fa. Note however, that at any higher order 
2k, k = 3,4... the weakest anisotropy of n.n. hopping will cause the same effects, in the 
effective Hamiltonians as well as in the corresponding phase diagrams. This implies that in 
the quantum systems described by H , any anisotropy of n.n. hopping orients the dimeric 
and axial-stripe phases in the direction of a weaker hopping. 

5 Summary 

We have considered two systems of correlated quantum particles, described by extended 
Falicov-Kimball Hamiltonians, with the hopping particles being either spinless fermions or 
hardcore bosons. The both system have been studied in the regime of strong coupling and 
half-filling, where the stability of some charge-stripe phases can be proved 0. Two main 
conclusions have been drawn. Firstly, any anisotropy of nearest-neighbor hopping orients the 
dimeric and axial-stripe phases in the direction of a weaker hopping. Secondly, even a weak 
anisotropy of hopping reveals a tendency of fermionic phase diagrams to become similar to 
the bosonic ones. 
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Appendix A 

Here we present some details pertaining to the m-potential method, which is used in this 
paper to construct ground-state phase diagrams. 

In comparison to the isotropic case, the presence of the hopping anisotropy lowers the 
symmetry of the Hamiltonians. Therefore, the zero-potentials K^? and the set of T-plaquette 
configurations have to be modified suitably. The zero-potentials, satisfying 

$>^ = 0, (28) 

T 

can be chosen in the form: 

9 

K^ = J2^4\ (29) 

i=l 

where the coefficients aj have to be determined in the process of constructing a phase di- 
agram, and the potentials k!^\ invariant with respect to the spatial symmetries of H and 
fulfilling condition (J2EJ), read: 

fop = Si + s 3 + s 7 + s g - 4s 5 , 
= S2 + Sg — 2^5, 
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(3) 


= S4 + Sq — 2s 5 , 




(4) 

rCrj-i 


= S1S2 + + S7S8 + SgSg — 2S4S5 






= + S^Sq + S4S7 + SqSq — 2S2S5 


— 2S5S8, 


(6) 


= + S3S5 + S 5 S 9 + S5S7 - S 2 S 4 - 


- s 4 s 8 - SgS 6 - s 2 s Gl 


(7) 


= S1S3 + S7S9 - 2s 4 s 6 , 






= S1S7 + S3S9 - 2s 2 s 8 , 






= S1S2S4 + S 6 S 8 Sg + S 2 S 3 S e + S 4 S 7 S 8 


- S2S4S5 - s 5 s 6 s 8 - s 2 s 5 s 6 - s 4 s 5 s 8 



In the above expressions the sites of a T-plaquette have been labeled 1, ... ,9, from left to 
right, starting at the bottom left corner and ending in the upper right one. In comparison 
with the isotropic case, the set of T-plaquette configurations is considerably larger and 
consists of 168 configurations (see Fig. [HJ). 

It turns up, that the coefficients «j can be chosen as affine functions of the parameters 
that enter linearly into the truncated effective Hamiltonian, i.e. energy parameters and 
the anisotropy parameter. Then, the energies of T-plaquette configurations become affine 
functions of these parameters. Consequently, the phase domains, being the solutions of a 
finite number of weak inequalities between the T-plaquette energies, are polyhedral convex 
sets in the space of energy and anisotropy parameters. Any point of a polyhedron is a convex 
combination of a finite number of points and directions (generating points and directions) 
|21] • The coefficients aij determined at the generating points and half-lines being boundaries 
of unbounded domains, enable one to find their values at any other point by taking suitable 
convex combinations. This is the content of the tables presented in Appendix C. 

Appendix B 

Here we provide some details referring to the presented phase diagrams: the sets of admissible 
T-plaquette configurations, used in constructing phase diagrams, equations of line boundaries 
between the phase domains, and coordinates of the crossing points of the line boundaries. 
The symbols like <Si|c>2 stand for the line boundary between the phases Si and S 2 , etc. 
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Figure 14: All the T-plaquette configurations in the case of anisotropic interactions, up to rotations by n and reflections 
in lattice lines parallel to the axes. The configurations that can be obtained from the displayed ones by the hole-particle 
transformation are not shown, only the numbers assigned to them are given in the brackets. 
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Fig UX6 15 '. All the admissible T-plaquette configurations, up to rotations by 7T, that are used in constructing the fourth-order 
phase diagram in a neighborhood of point a of Fig. HI i.e. the elements of 




1 2 3 4 5 6 

Fig UX6 16 i All the admissible T-plaquette configurations that are used in constructing the fourth-order phase diagram, up 
to rotations by n, in a neighborhood of point b of Fig. EI i- e - the elements of 5y (the vertical n.n. pairs that are occupied by 
one ion are forbidden). 
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Figure 17: All the admissible T-plaquette configurations, up to rotations by it and reflections in lattice lines parallel to the 
axes, that are used in constructing the fourth-order phase diagram in a neighborhood of point c~ of Fig.|3 i.e. the elements 
of 5y (the horizontal n.n. pairs that are occupied by two ions are forbidden). 
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Table 1: Domain boundaries of the phase diagrams for fermions, shown in Fig. ^ 
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u = fi2 + l 


•Ssegl'Sdd 


e = 2u - 2(3 2 + 3 




■Sdd\Scb 


e = -2w + 2/3 2 + 21 




Sseg\S%l 


e = u; + ± 


•Sdll^cb 


e = -u + 2/3 2 + f 




cu = 2/? 2 - | 




Sdd\Sdi 


a; - T 






£ = 2/? 2 - 6 




cv 1 cv 
C> vl\ v2 


e = w-| 


Cv 1 cv 
^l^dl 


e = -w + 2/3 2 + § 


Cv 1 cv 




co = (3 2 + | 


Points 


< 02 < 7 


7 < 02 < 10 


a 


(A + §,12) 


(/3 2 + |,/3 2 + 5) 


b 


(2/3 2 - |,2/3 2 -2) 


(/% + §,& + !) 


c 


(f,2ft-2) 




d 


(2/3 2 -|,2/3 2 -6) 




e 


(f,2/? 2 -6) 
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Table 2: Domain boundaries of the phase diagrams for hardcore bosons, shown in Fig. ^ 



Lines 


< 2 2 


2 < 02 < 5 


Sseg | Set 


w = 02 + | 


Sseg | Spot 


e = 2a; - 20 2 + 4 




Spcb | Scb 


e = -2w + 20 2 + 14 




c 1 cv 
°seg\O vl 


e = u; + ! 


Sdi\S c b 


e = -w + 20 2 + f 


S^i\Spcb 


w = 20 2 + | 




Spcb\Sdl 


w = | 






e = 20 2 + 1 




CV 1 cv 
<\ll<-\;2 


e = u; + i 




e = -w + 20 2 + § 






W = 02 + | 


Points 


< 2 «S 2 


2 < 2 < 5 


a 


(A + |,9) 


(02 + 1,^2+7) 


b 


(20 2 + ±,20 2 + 5) 


(02 + 1,02+3) 


c 


(§,202 + 5) 




d 


(20 2 + |,20 2 + l) 




e 


(§,20, + 1) 
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Table 3: Domain boundaries of the phase diagrams shown in Fi, 



Lines 


< 02 < 1 


K 02 sS 2 


2 < /3 2 < 3 


3 < 02 sS 6 


S+\S- 


5 = 


Sdd\S- 


X — W i 02 3 
2^2 4 




Sdd\Sl 


X — U , 02 11 
6 1 6 12 




Sdd\S2 


* = -f 




Sdd\S3 


£ _ W 02 21 
2 2 4 




S-\S X 


5 = -u) + /3 2 - i 


Sl\S 2 


<5 = -w + 2 + 2 


5 2 |<S3 


(5 = -tv + 2 + 6 




<5 = + /3 2 + f- 






c a; 1 
5 = ~4 ~ 8 






° 4^2 8 




Sd 2 \S- 




r _ 3w , /?2 3 
° — 4^2 8 






x _ w 3/3 2 , 1 
° — 4 2 8 






r _ w ?2 , 27 
° — 4 2 8 


\$dd 




W = 202 - | 




S dl\ S dd 




, i — 23 




S dl \ S cb 




uj = 20 2 + ±§ 


CU 1 ov 
^Pul 




° 4 8 


S^ 3 \S d d 




° 2^2 4 




W- 




x w i 1 
5 = -2 + 4 


CD 1 CD 




° 4 2^8 






c _ 3W , 02 1 13 

° — 4 ^ 2 "i" 8 


S d4,\ S dd 




° 2 2 








5 = -#2 + 5 


CD 1 CD 




c _ 3u , i 37 
° ~~ 4 2 8 






X — w , §2 29 
° — 2 ^ 2 4 




s d:j\ s 3 




5 = | - 2/3 2 - | 


CD 1 CD 

^dil^di 




c _ 3u , /3 2 , 57 
° — 4 1" 2 8 


CD 1 CD 

C3 v2\ C:, v3 




5 = - 2 + 4 


CD 1 CD 

C 'v2\ C 'd4 




<5 = -/3 2 + 2 


CD 1 CD 




5 = | - 202 + ± 


CD 1 CD 




W = 202 + ^ 
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r OlIICS 


< /3 2 < 1 


1< /3 2 < 2 


2 < ^ 3 


3 < 2 ^ 6 


a 


(A -|,o) 


H>°) 


D 


(A + 1,-1) 


(2#2 - |, ~t + 1) 


(2/32 - |, t " 1) 






(A +1,-1) 


(23 ft , i\ 

V 2 ' 2 2/ 


(23 ft 3\ 
V 2 ' 2 2/ 


(3 _l) 

V2' 2/ 


d 


(& + ¥,-§) 


(2/3 2 + f,-/3 2 + l) 


(i.-i) 


(2/3 2 + f , -/3 2 + 3) 


e 


(/3 2 + f,0) 


(3 ft _ 3\ 
V2' 2 2/ 


(4/3 2 -f,-f + 1) 


(2/3 2 -|,-/3 2 + 2) 


f 




(4^-1,-f-l) 


(2/3 2 -|,-/3 2 + |) 


(2/3 2 + |,-/3 2 + 2) 


g 




(2/3 2 -|,-/3 2 ) 


(2/3 2 + f ,-/3 2 + |) 


(2/3 2 -|,-/3 2 + |) 


h 




(A + I> -1) 


(11 ft 5\ 
V 2 ' 2 2/ 


(2/3 2 + f , -/3 2 + |) 


i 




&% + ¥.-§) 


(19 ft 5\ 
I 2 ' 2 2/ 


(2/3 2 + f,-/3 2 + l) 


j 




(2/3 2 + f,-/3 2 + l) 


(2/3 2 + |,-/3 2 + |) 


k 




(2/3 2 + |,-/3 2 + i) 


(2/3 2 + -/3 2 + |) 


1 




(2/3 2 + ^,-/3 2 + i) 


(2/3 2 -i,-/3 2 ) 


m 




(2/3 2 -|,-/3 2 ) 
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Table 4: Domain boundaries 



of the phase diagrams shown in Fig. El 



Lines 


< /3 2 ^ 3 




5 = 




A- _y _ £ 
4 8 




5 = - + - 


cv 1 cv 


(5 — -- - i 

2 4 




A — _^ _ 1 
2 4 


cv 1 cf 


r w 3ft 23 
— 4 2 8 


CV 1 CD 

D d 2 l D d3 


X _ 3a; , ft i 5 
u — 4 "•" 2 "i 8 




£ _ 3uj , ft 11 
~ 4 "i" 2 8 


5 d 2 l 5 l 


r a; 3ft 39 
— 4 2 8 


CD 1 CD 

D D2l D d3 


(5 = f - 2/3 2 - f 




(5 = f - 2/3 2 - f 


CD 1 CD 

D d3l D dl 


: 3oj ■ ft ■ 25 

— 4 2 8 




I w ft 5 

u 4 2 8 


<Wi 


W = 2/3 2 + ^ 




cj = 2/3 2 + f 




<5 = -w + /3 2 - | 


5i|5 3 


5 = -u + p 2 + 2 


^3 5cb 


6 = -u + /3 2 + f 



Points 


< /3 2 ^ 3 


a 


(-l.o) 


b 


(-|,o) 


c 


V 2 ' 2 / 


d 


(2/3 2 + 4±,-/3 2 -l) 


e 


(2/3 2 + l,-^-2) 


f 


(2/3 2 + |,-^-l) 


g 


(2^ + f,-^2-|) 


h 


(2/3 2 + f ,-/3 2 -3) 


i 


(2/3 2 + ^,-/3 2 -7/2) 


j 


(2^ + 1,-/32-4) 
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Table 5: Domain boundaries of the phase diagrams shown in Fig. |H] and Fig. 



Lines 


fermions 


bosons 


^seg I £cb 


i , — 9 


w = | 




9 


w = | 




e = w + i 


e = w + | 




e = - W + ¥ 


cv 1 cv 


e = w - | 


e = w + i 


CV 1 ci> 


e = ~ L0 + T 



Points 


fermions 


bosons 


a 


(1,5) 


(1,7) 


b 


(1.1) 


(1,3) 



Table 6: Domain boundaries of the phase diagrams shown in Fig. IHland Fig. El 



Lines 


fermions 


bosons 








s-\s+ 


5 = 


= 




Points 


fermions 


bosons 


s v vl \s- 


6 = -- - - 

u 4 8 


u 4 8 




a 


H,o) 


(-l.o) 


Cv 1 cv 


x - a - 1 

4 8 


5 — - + - 

4^8 




b 


(i.o) 


(-io) 




(5 = — - + \ 

u 2^4 


u 2 4 




c 


V2> 2/ 


f_5 _n 

V 2' 2/ 


Cv 1 cv 


2^4 


u 2 4 









29 



Table 7: Domain boundaries of the phase diagrams shown in Fig. [HI] and Fig. EI 



Lines 


fermions 


bosons 


d d3 


S- 


6=--+- 
2^4 


w 2 4 


cv 
D d3 


cv 
D d2 


u 2^4 


5 = -- - - 

° 2 4 


cv 
^v3 


CD 

D d2 


S = OL + 11 
° 4^8 


X _ 0) 23 
~~ 4 8 


cv 
D d2 


S- 


s: 3a; 3 
~~ 4 8 


r 3w 11 
— 4 8 


cv 
D d2 


CD 


r 3uj 1 13 

u ~ 4 8 




cv 
^d2 


CD 

D d3 




A = _3w , 5 
u 4^8 


cv 
D d2 


«5i 


5 = - + - 

u 4^8 


r w 39 
— 4 8 


cv 


CD 

D d2 


5 = 2 




cv 


5 2 


8 = \ 






CD 

D d3 


s. 3iii 1 37 
~ 4 "r" 8 




CD 


CD 

D d2 


5 — - + - 

2 -1- 4 


A — w _ I! 
"24 


CD 

D d3 


S 3 


u 2 4 


r £J 25 
— 2 4 


cv 
D d3 


CD 


r 3u 1 57 
— 4 ~r 8 


r 3u 1 25 
— 4 " r 8 


CD 
<->dl 


S 3 


4^8 


A — w _ 5 
4 8 


PI) 


CD 

<-\;2 




cv 






S-\S! 


5 = -w - i 


S!\S 2 


5 = -lo + 2 




Si\S 3 




5 = -u + 2 




^ = -co + 6 




S3 \S c b 


<5 = - w +f 


J = - u + M 



Points 


fermions 


bosons 


a 


(¥.3) 




b 


H,2) 


(1,-2) 


c 




(3 _5\ 
V2> 2/ 


d 


f_5 3\ 
V 2' 2/ 


(¥,-1) 


e 


(f,l) 


(* -3) 


f 


(f>l) 


en _i) 

V 2 ' 2/ 


g 


f3 1\ 
V2> 2/ 


(1,-4) 


h 


fll I") 

V 2 ' 2/ 




i 


H.°) 
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Appendix C 

Below, in a series of tables, we provide the coefficients {a,}, i = 1, . . . , 9, of the zero-potentials 
for the phase diagrams presented in this paper. The coefficients that are missing in a table 
are equal to zero. As in Appendix B, the symbol <Si|<S2 denotes the boundary between the 
phases S\ and S2 , etc. 

Table 8: Zero-potentials coefficients for the phase diagram of fermions, shown in Fig. Q for 
fa = 0. 





•Sseg\<Scb <S seg \S v i S c b\Sdl 
















a 3 





04 


CO 1 3 

96 ~r 64 


tu 1 7 
96 64 


w 1 
96 64 


a 5 


LU 1 3 

96 ~r 64 


w 1 7 
96 64 


w 1 
96 64 


a 6 





1 

8 




1 

48 


1 

12 


a 8 


1 

48 


1 

12 



Table 9: Zero-potentials coefficients for the phase diagram of fermions, shown in Fig. ^ for 
ft = 7/2- 





C 1 C CI C C 1 QV cv 1 cv cv 1 c» 
OseglOcf, £> seg \£> vl O cb \£> dl O v2 \O vl O v2 \O dl 


«1 





«2 





«3 





«4 


u> I 3 
96 64 


«5 


u , 23 
96 " r 192 


«6 





OL7 


1 

48 


«8 


1 

48 
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Table 10: Zero-potentials coefficients for the phase diagram of fermions, shown in Fig. ^ for 

ft = r. 



S n, \ 


C 1 C CI cv C 1 cv 
<~>seg\<~>cb c> seg\ c > v l °cft| < - , rfi 


CI) 1 CD C" 1 C 

<\2l°ul D i>2l°dl 


Ot\ 





Oi2 





«3 





«4 


U) 1 

96 64 


w 1 3 
96 64 


a 5 


w , 49 
96 192 


lo 1 37 
96 " r 192 


«6 


1 

8 





OL-J 


1 

12 


1 

48 


a s 


1 

12 


1 

48 



Table 11: Zero-potentials coefficients for the phase diagram of fermions, shown in Fig. ^ for 
P2 = 10. 



{"J 


c 1 c CI cv C 1 C w 
^seg^cb c> seg\ i -' v l i - > cb\ c ' c [i 


cv I cv cv | cv 
<-V2l<-V <\>2l°dl 







«2 





«3 





«4 


1 

6 


u) , 3 
96 64 


«5 


1 

6 


w , 49 
96 ~r 192 


«6 


1 

8 





«7 


1 

12 


1 

48 


«8 


1 

12 


1 

48 
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Table 12: Zero-potentials coefficients for the phase diagram of hardcore bosons, shown in 
Fig. □ for 2 = 0. 



{<*%! 


•Sseg\<Scb S se g\S v x S c b\Sdl 


<Si;2 |<Sul 




Ot\ 





Oil 





a 3 





a.4 


LU i 5 

96 ' 192 


w i 17 
96 " r 192 


u> 7 
96 192 


a 5 


LU i 5 

96 " r 192 


w , 17 
96 " r 192 


ui 7 
96 192 


«6 





1 

8 


aj 


1 

48 


1 

12 


a 8 


1 

48 


1 

12 



Table 13: Zero-potentials coefficients for the phase diagram of hardcore bosons, shown in 
Fig. [U for (3 2 = 1. 



{"J 


Sseg | <5cb 


c 1 cv c 1 c" Qv 1 cti cv 1 cv 
&seg\<-> v l ^>cb\O dl O^l^l <-\,2l D dl 







«2 





«3 





«4 


1 

24 


lu i 5 
96 ~r 192 


«5 


1 

24 


oi i 3 
96 " r 64 


«6 





«7 


1 

48 


«8 


1 

48 
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Table 14: Zero-potentials coefficients for the phase diagram of hardcore bosons, shown in 
Fig. □ for (3 2 = 2. 





Sseg | S c b 


C 1 cv 
°seg\O v i 


S c b\S v dl 


cv 1 cv 


CV 1 cv 
D i>2Pdl 







a 2 





«3 







1 

12 


u> 1 
96 192 


U! 1 1 

96 24 


«5 


1 

12 


ui , 19 
96 ' 192 


a 6 


1 

8 


1 

32 


1 

32 


Oi-J 


1 

12 


1 

48 


1 

192 




1 

12 


1 

48 


5 

96 



Table 15: Zero-potentials coefficients for the phase diagram of hardcore bosons, shown in 
Fig. d for = 5. 





C 1 C CI C 1 CV 
t-'seg^cb c> seg\<~> v l c, cb\ c ' c li 


CV 1 CV CV 1 ct> 
^l^vl <-\>2l D cZl 


ai 





«2 





a 3 







11 

96 


US i 1 

96 24 


«5 


11 

96 


ui i 19 
96 ~r 192 


«6 


1 

8 


1 

32 




1 

12 


1 

192 


a 8 


1 

12 


5 

96 
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Table 16: Zero-potentials coefficients for the phase diagram shown in Fig. |2l for /3 2 





S + \S- 


S-\S! 


5i|5 2 


<?2 ^3 


S3 \$cb 


5 = 0, uj > f 







UJ | 3 

16 ' 32 


uj 3 
16 32 


u 9 
16 32 


uj 21 
16 32 










UJ 1 1 

8 16 


3 

8 16 


uj 3 
8 4 


ui 21 
8 16 





a 3 





U> 1 1 

8 ~r 16 


cj 3 
8 16 


cj 3 
8 4 


w 21 
8 16 





«4 


UJ 1 3 

96 ~r 64 


u 1 1 
96 ' 16 




cj i 3 
96 " r 64 


a 5 


w i 3 
96 ~r 64 


w . 1 
96 "r" 16 




i 3 
96 64 


«6 





1 

16 


1 

32 





OL7 


1 

48 


1 

96 


1 

48 


«8 


1 

48 


1 

96 


1 

48 



Table 17: Zero-potentials coefficients for the phase diagram shown in Fig. for fa 







S-\S! 


«5i|5 2 


<?2 £3 


S-i\Sdb 


5 = 0, u > f 


o>\ 





U> 1 1 

16 " r 32 


uj 3 
16 32 


uj 11 
16 32 


uj 23 
16 32 










UJ 1 

8 16 


W 11 

8 32 


uj 41 
8 48 


uj 23 
8 16 





«3 





UJ 1 

8 16 


uj 13 
8 32 


uj 43 
8 48 


uj 23 
8 16 





«4 


uj i 3 
96 " r 64 


uj _,_ 9 
96 ~r 128 


uj 1 13 
96 ~r 192 


uj 1 3 
96 ~r 64 


a 5 


w 1 13 
96 1 192 


uj | 29 
96 1 384 


uj 1 5 
96 ~r 64 


uj 1 13 
96 ~r 192 


a 6 





1 

32 





07 


1 

48 


1 

96 


1 
192 


1 

48 


a 8 


1 

48 


1 

96 


5 
192 


1 

48 



35 



Table 18: Zero-potentials coefficients for the phase diagram shown in Fig. El for /3 2 = 3/2. 



lines & points 


a\ 


«2 


a 3 


«4 


«5 


a 6 


a-j 


«8 


a 9 













iii 1 3 
96 64 


ui 1 5 
96 " r 64 


l 

64 


7 
192 


1 

48 





S-\S! 


16 


u 1 

8 8 


uj 3 
8 32 


uj 1 15 
96 256 


u) | 5 
96 64 


1 

256 


17 

384 


1 
192 





5i|5 2 


O) 49 
16 256 


u> 31 
8 64 


uj 13 
8 32 


ui | 7 
96 1 128 


w , 27 
96 1 256 


7 
256 


23 
384 


23 
384 


19 
256 


<$2 M>3 


uj 105 
16 256 


uj 115 
8 128 


uj 29 
8 32 


uj _i_ 7 
96 ~ l ~ 128 


u> , 43 
96 ~ l ~ 512 


25 
512 


157 
3072 


113 
3072 





Ss\Scb 


UJ 3 

16 4 


ui 3 
8 2 


uj 3 
8 2 


u) i 3 
96 64 


i 5 
96 " r 64 





1 

48 


1 

48 





5 = 0, uj > f 











uj 1 3 
96 64 


i 5 
96 " r 64 


1 

64 


1 

48 


7 
192 





(|,o) 











1 

24 


7 

96 





1 

48 


1 

48 





(¥>o) 











7 

96 


1 

24 





1 

48 


1 

48 





b 


1 

32 





1 

16 


1 

24 


7 

96 





1 

48 


1 

48 





c 


1 

32 


1 

16 





7 

96 


1 

24 





1 

48 


1 

48 





e 


3 

32 


1 

16 


1 

8 


1 

32 


1 

16 





1 

48 


1 

48 





f 


7 

64 


1 

8 


3 

16 


1 

96 


11 
192 


1 

64 


5 

96 


5 

96 


3 

64 



Table 19: Zero-potentials coefficients for the phase diagram shown in Fig. |21 for (3 2 — 2. 



lines & points 




a 2 


a 3 


«4 


a 5 


o 6 


OL-J 


a s 


«9 


S + \S- 











D I 3 
96 64 


uj , 17 
96 192 


1 

32 


5 

96 


1 

48 





S-\Si 


UJ 1 

16 32 


uj 3 
8 16 


uj 3 
8 16 


uj | 3 
96 64 


w , 17 
96 192 





1 

48 


1 

48 





5i|5 2 


UJ 1 

16 4 


uj 9 
8 16 


uj 7 
8 16 


uj , 3 
96 " r 64 


0) , 23 
96 192 


1 

32 


1 

12 


1 

12 


3 

32 


S2 <?3 


uj 15 
16 32 


uj 15 
8 16 


uj 15 
8 16 


| 3 
96 ~r 64 


w , 17 
96 192 


1 

16 


1 

12 


1 

24 





<?3 5cfc 


ui 25 
16 32 


uj 25 
8 16 


uj 25 
8 16 


D 1 3 
96 ~r 64 


w , 17 
96 ~r 192 





1 

48 


1 

48 





5 = 0, UJ > f 











1 3 
96 ~r 64 


w , 17 
96 192 


1 

32 


1 

48 


5 

96 





(l,o) 











1 

32 


7 

96 





1 

48 


1 

48 





(f,o) 











7 

96 


1 

32 





1 

48 


1 

48 





b(e) 


1 

16 





1 

8 


1 

32 


7 

96 





1 

48 


1 

48 





c 


1 

16 


1 

8 





7 

96 


1 

32 





1 

48 


1 

48 
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Table 20: Zero-potentials coefficients for the phase diagram shown in Fig. El for 02 = 5/2. 



lines & points 


OL\ 


«2 


a 3 


«4 


«5 


a 6 


017 


a 8 


a 9 


S+\S- 











uj 1 3 
96 64 


uj 1 19 
96 192 


1 

32 


5 

96 


1 

48 





S-\St 


UJ 

16 16 


UJ 

8 4 


UJ J. 

8 4 


U 1 O 

96 "•" 64 


/ . in 
UJ 1 ly 

96 192 





1 
i 

48 


.1 

48 





Si\S 2 


oj y 
16 32 


UJ D 

8 8 


UJ 1 
8 2 


UJ \ 

96 "•" 64 


a; 1 zo 
96 192 


1 

32 


1 

12 


1 

12 


3 

32 


<?2 ^3 


w 33 
16 64 


uj 31 
8 32 


- — 1 
8 1 


uj 1 3 
96 ~r 64 


w 1 35 
96 384 


9 
128 


1 

12 


C 


192 





<?3 \$cb 


/1 1 Q 

CJ lo 

16 16 


/1 1 Q 

8 8 


< . 1 Q 

UJ lo 

8 8 


UJ 1 O 

96 ^ 64 


cj _i_ iy 
96 192 





1 

_ 48 


— 48 





5 = 0, to > f 











1 

96 "•" 64 


uj i iy 
96 192 


32 


-I 

_ 48 


C 



96 





(i.o) 











1 

48 


7 

96 





1 

48 


1 

48 





(f>o) 











7 

96 


1 

48 





1 

48 


1 

48 





b 


1 

16 


1 

4 


1 

4 


5 

96 





3 

32 


1 

12 


1 

48 





c 


1 

32 


1 

16 





7 

96 


1 

48 





1 

48 


1 

48 





d 


1 

16 





1 

8 


1 

32 


1 

12 





1 

48 


1 

48 





e 


1 

32 


1 

8 


3 

16 


1 

96 


1 

16 


1 

16 


5 

96 


5 

96 





f 


1 

8 


1 

8 


1 

8 


7 
192 


7 

96 


1 

64 


1 

96 


1 

48 





g 


3 

32 


3 

16 


1 

16 


1 

12 


3 

64 


1 

64 


1 

48 


1 

48 





h 


1 

16 


1 

8 


1 

4 


1 

96 


7 

96 


1 

16 


1 

12 


1 

12 


1 

16 


i 


1 

16 


1 
4 


1 
4 


5 

96 





3 

32 


1 

12 


1 

48 
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Table 21: Zero-potentials coefficients for the phase diagram shown in Fig. |21 for /?2 = 3. 



lines & points 




a 2 


«3 


04 


a 5 


a 6 


«7 


a 8 


S+\S- 











w 1 3 
96 64 


w i 7 
96 64 


i 

32 


5 

96 


l 

48 


S-\St 


UJ 3 

16 32 


u 5 
8 16 


w 5 
8 16 


u 1 3 
96 64 


ti> 1 7 
96 64 





1 

48 


1 

48 


Si\S 2 


ui 11 
16 32 


uj 7 
8 16 


UJ 11 
8 16 


£J | 3 

96 64 


ct> 1 3 
96 64 


1 

16 


1 

12 


1 

12 


s 2 \s 3 


uj 9 
16 16 


- - 1 

8 


uj 17 
8 16 


O) i 3 
96 64 


k> i 3 
96 32 


5 

64 


1 

12 


1 

96 


<?3 \S c b 


oj 27 
16 32 


w 27 
8 16 


uj 27 
8 16 


cj 1 3 
96 64 


u> , 7 
96 ~r 64 





1 

48 


1 

48 













ixJ 1 3 
96 ' 64 


w i 7 
96 ~r 64 


1 

32 


1 

48 


5 

96 


(1,0) (b) 











1 

96 


7 

96 





1 

48 


1 

48 


(f ,0) (c) 











7 

96 


1 

96 





1 

48 


1 

48 


d 


1 

16 





1 

g 


1 

32 


3 
32 





1 

48 


1 

48 


e(h) 





1 

4 


1 

4 


1 

96 


5 

96 


1 

8 


1 

12 


1 

12 


i 





1 

4 


1 

4 


5 

96 


1 

96 


1 

8 


1 

12 


1 

12 


f 


5 

32 


7 

32 


3 

16 


35 
768 


7 

96 


3 
256 


5 
384 


1 

48 


g 


1 

8 


1 

4 


1 
8 


3 

32 


3 

64 


1 

64 


1 

48 


1 

48 



Table 22: Zero-potentials coefficients for the phase diagram shown in Fig. |21 for f3 2 — 6. 



lines & points 


ati 


a 2 


a 3 


«4 


a 5 


«6 


aj 


a s 


S+\S- 











uj 1 3 
96 ~r 64 


UJ , 11 
96 64 


1 

32 


5 

96 


1 

48 


S-\Si 


uj 9 
16 32 


ui 11 
8 16 


UJ 11 

8 16 


uj I 3 
96 64 


u , 11 
96 64 





1 

48 


1 

48 


5i|5 2 


u> 19 
16 32 


UJ 11 
8 16 


ui 15 
8 16 


uj 1 3 
96 "•" 64 


O) 1 7 
96 " r 64 


1 

8 


1 

12 


1 

12 


S 2 \Ss 


ui 25 
16 32 


ui 21 
8 16 


ui 23 
8 16 


uj I 3 
96 64 


£J , 9 
96 ~r 64 


3 

32 


1 

12 


1 

48 


S3 \S c b 


uj 33 
16 32 


ui 33 
8 16 


uj 33 
8 16 


uj 1 3 
96 64 


UJ , 11 
96 ~r 64 





1 

48 


1 

48 













uj 1 3 
96 ~r 64 


UJ , 11 
96 64 


1 

32 


1 

48 


5 

96 


b 











1 

96 


13 
96 


1 

32 


5 

96 


1 

48 


(f,o) 











13 
96 


1 

96 


1 

32 


1 

48 


5 

96 


c 


1 

16 





1 

8 


1 

32 


5 

32 





1 

48 


1 

48 


d 


3 

16 


3 
8 


3 
8 


13 
96 


5 
192 


1 

64 


1 

48 


1 

48 


e 


7 

32 


5 
8 


9 
16 


11 
192 


5 

96 


7 

64 


1 

12 


1 

12 


f 


3 

16 


5 
8 


5 
8 


11 
96 


1 

48 


3 

32 


1 

12 


1 

48 


g 


11 

32 


5 
8 


9 
16 


1 
192 


7 

96 


1 

64 


1 

48 


1 

48 


h 


5 

16 


5 
8 


1 

2 


5 

32 


3 

64 


1 

64 


1 

48 


1 

48 



38 



Table 23: Zero-potentials coefficients for the phase diagram shown in Fig. 01 for f3 2 — 0. 



lines & points 




a 2 


a 3 


«4 


a 5 


«6 


OL-j 


a 8 


S+\S- 











u> i 5 
96 " r 192 


uj , 5 
96 192 





1 

48 


l 

48 


S-\St 


ui i 3 
16 " r 32 


UJ _i_ 1 

8 16 


UJ i 1 

8 "r" 16 


oj i 5 
96 " r 192 


£J i 5 
96 192 





1 

48 


1 

48 


St\S 3 


uj 9 
16 32 


ui 1 
8 16 


UJ 1 

8 16 


ui 7 
96 192 


w 7 
96 192 


1 

8 


1 

12 


1 

12 


<?3 \Scb 


a; 13 
16 32 


ui 13 
8 16 


uj 13 
8 16 


cj 1 5 
96 192 


uj i 5 
96 192 





1 

48 


1 

48 


5 = 0, u > f 











uj , 5 
96 " r 192 


oj i 5 
96 ' 192 





1 

48 


1 

48 


b 











3 

32 


3 

32 


1 

8 


1 

12 


1 

12 













3 

32 


3 
32 


1 

8 


1 

12 


1 

12 


c 


3 

32 


1 

16 


1 

16 


1 

12 


1 

12 


1 

32 


1 

96 


1 

96 


d 


1 

16 


1 

8 


1 

8 


3 

32 


3 

32 


1 

8 


1 

12 


1 

12 


e 


1 

16 


3 
8 


3 
8 


1 

96 


1 

96 


1 

8 


1 

12 


1 

12 


f 


3 

32 


9 
16 


9 
16 


7 

96 


7 

96 


3 
32 


1 

12 


1 

12 


g 


5 

32 


5 

16 


5 

16 


1 

12 


1 

12 


1 

32 


1 

48 


1 

48 



Table 24: Zero-potentials coefficients for the phase diagram shown in Fig. 01 for (3 2 — 3. 



lines &; points 


oil 


a 2 


«3 


«4 


a 5 


«6 


0/-J 




S+\S- 











ui , 5 
96 ~r 192 


ui i 17 
96 ~r 192 


1 

32 


5 

96 


l 

48 


S-\Sx 


uj 3 
16 32 


uj 5 
8 16 


uj 5 
8 16 


ui , 5 
96 192 


ui , 17 
96 ~r 192 





1 

48 


1 

48 


Sx\S 3 


uj 15 
16 32 


uj 7 
8 16 


uj 7 
8 16 


uj 7 
96 192 


ui 1 5 
96 192 


1 

8 


1 

12 


1 

12 


<?3 \$cb 


uj 19 
16 32 


u> 19 
8 16 


ui 19 
8 16 


u> 1 5 
96 ~r 192 


ui , 17 
96 ~r 192 





1 

48 


1 

48 


5 = 0, U^f 











ui i 5 
96 192 


ui i 17 
96 192 


1 

32 


1 

48 


5 

96 


b 











1 

32 


3 

32 


1 

32 


5 

96 


1 

48 


(¥,o) 











3 

32 


1 

32 


1 

32 


1 

48 


5 

96 


c 


3 
32 


1 

16 


1 

16 


13 
192 


25 
192 


1 

64 


1 

96 


1 

48 


d 


1 

4 


1 

2 


1 

2 


3 

32 


3 

64 


1 

64 


1 

48 


1 

48 


e 


1 

4 


3 
4 


3 
4 


7 

96 


1 

96 


1 

8 


1 

12 


1 

12 


f 


13 
32 


23 
32 


11 
16 


5 
384 


1 

96 


1 

128 


1 

48 


1 

48 


g 


11 

32 


11 
16 


11 
16 


7 

48 


13 
192 


1 

64 


1 

48 


1 

48 



39 



Table 25: Zero-potentials coefficients for the phase diagram shown in Fig. |HJ 



lines 


Oil 









U) 1 1 

16 32 


cv 1 cv 


w 7 
16 32 


8 = 0, uj^ I 






Table 26: Zero-potentials coefficients for the phase diagram shown in Fig. El 



lines 


ati 


S+\S- 





S-\S& 


oj , 9 
16 ~r 32 


Cv 1 cv 
D d3I D d2 


LU _|_ 1 

16 32 


8 = 0, u > -\ 






Table 27: Zero-potentials coefficients for the phase diagram shown in Fig. El 



lines & points 




«2 


«3 


Ct4 


a 5 


a 6 


«7 


a 8 


s-\s v v3 








LU I 1 

8 ~r 16 


a; | 1 
24 16 


LU 1 3 

96 64 


1 

32 


1 

48 


l 

48 


S-\Si 





1 

8 


UJ 1 1 
4 T 4 


7oj 3 
96 64 


io 1 3 
96 t 64 





1 

48 


1 

48 


Si\S 2 


7^ 17 
96 64 


7u> 1 
48 32 


5oj 3 
48 32 





CJ 1 

96 64 


1 

8 


1 

12 


1 

12 


52^3 


7oj 29 
96 64 


7lu 21 
48 32 


5a; 19 
48 32 





a; 1 1 
96 ~r 64 


3 

32 


1 

12 


1 

48 


^3 \S c b 








ui 21 
4 8 


7w , 45 
96 "r" 64 


LU 1 3 

96 t 64 





1 

48 


1 

48 










5 
4 


5 5 
16 96 


5 

96 





1 

48 


1 

48 


cv 1 ci) 
C, v2\ C> dl 








4 


S 1 
16 96 


5 

192 


1 

64 


1 

48 


1 

48 


cv 1 cv 
C> v3\ v2 


1 

8 





LU 3 

8 16 


LU | 3 

24 t 32 


LU | 3 

96 ' 64 


1 
8 


1 

12 


1 

12 


C 


1 

8 





1 

4 


7 

96 


1 

48 


3 

32 


1 

12 


1 

48 


e 








3 

8 


1 

32 


3 

64 


1 

64 


1 

48 


1 

48 



40 



Table 28: Zero-potentials coefficients for the phase diagram shown in Fig. ITT 



lines &; points 




«2 


«3 


04 


a 5 


a 6 


OL-j 


a 8 


s-\s v v3 





1 

4 


w i 9 
8 ~r 16 


LU 7 

24 48 


LU I 5 

96 " r 192 


i 

32 


1 

48 


1 

48 


S-\St 





1 

8 


Ul 1 1 

4 "i 4 


7a; 13 
96 192 


u i 5 
96 " r 192 





1 

48 


1 

48 


Si\S 3 


l 

4 





0) 1 
4 8 


7a; 1 
96 192 


LJ 7 

96 192 


1 

8 


1 

12 


1 

12 


Ss\Scb 








a; 13 
4 8 


7a; , 83 
96 " r 192 


u i 5 
96 " r 192 





1 

48 


1 

48 










4 


5 7 
16 96 


7 

96 





1 

48 


1 

48 


cv 1 cv 








<5 
4 


5 1 
16 32 


3 

64 


1 

64 


1 

48 


1 

48 


cv 1 cv 
<\>3l D i>2 


1 

8 





oj i 5 
8 ~r 16 


a; 5 
24 96 


a) i 5 
96 " r 192 


1 

8 


1 

12 


1 

12 


d 








5 
8 


23 
96 


13 
192 


1 

64 


1 

48 


1 

48 



41 



